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ABSTRACT 

To develop future Reusable Launch Vehicle concepts, we have conducted design optimization for a single element 
rocket injector, with overall goals of improving reliability and performance while reducing cost. Computational solutions 
based on the Navier-Stokes equations, finite rate chemistry, and the k-e turbulence closure are generated with design of 
experiment techniques, and the response surface method is employed as the optimization tool. The design considerations 
are guided by four design objectives motivated by the consideration in both performance and life, namely, the maximum 
temperature on the oxidizer post tip, the maximum temperature on the injector face, the adiabatic wall temperature, and the 
length of the combustion zone. Four design variables are selected, namely, H 2 flow angle, H 2 and O 2 flow areas with fixed 
flow rates, and O 2 post tip thickness. In addition to establishing optimum designs by varying emphasis on the individual 
objectives, better insight into the interplay between design variables and their impact on the design objectives is gained. 
The investigation indicates that improvement in performance or life comes at the cost of the other. Best compromise is 
obtained when improvements in both performance and life are given equal importance. 
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1. INTRODUCTION 

The Space Shuttle has been propelled into space for the last 20 years by the world’s most sophisticated reusable engine, 
the Space Shuttle Main Engine (SSME). However, due to limited knowledge base and stringent operating conditions, the 
engine’s key hot-section components, including combustion devices and turbomachinery, have encountered reliability 
issues. Consequently, frequent maintenance and periodic redesign have been required, resulting in high operation costs. The 
goal of NASA’s Next Generation Launch Technology (NGLT) program is to develop technologies to make the next 
generation launch systems safer, more affordable and more reliable. Design of robust combustion devices, namely, 
injectors, chambers and nozzles, is an important element in meeting these ambitious goals. 

Thrust chamber is an important part of a liquid propellant rocket engine. It takes in the propellants through the injector, 
atomizes, vaporizes, mixes and bums it in the combustion chamber and accelerates the produced gases, which are then 
ejected out to produce the required thrust. Hence, the performance and the environment influencing it are, to a large extent, 
governed by the design of the injector. The performance is determined by the rate and the extent to which the propellants 
mix and resultant combustion occurs. This mixing efficiency is a function of the details of the injector design. On the other 
hand, the location of the mixing and combustion is a key in determining the thrust chamber thermal environment. This 
environment directly affects the component life, including temperatures on the combustor wall, injector face and, for 
coaxial injectors, the oxidizer post tip. The difficulty encountered in designing injectors that perform well and have 
manageable environment is that, generally, the factors that promote performance also lead to increased heat fluxes to the 
solid surfaces in the combustor. 

Current injector design tools have been in use for 30 years or more^'^, and are largely empirical-based. The 
experimental databases, and thus the tools developed from them, are limited, in terms of design space, to specific element 
configurations that have been tested^. In terms of scope, the design tools typically focus on performance with the 
environment, being a secondary consideration. The limited amount of environmental information available from these tools 
is usually sketchy and not calculated as a function of the details of the injector design. It is very doubtful that application of 
these traditional design tools will result in future propulsion devices that meet the goals noted above. 

With rapid progress made in computational fluid dynamics (CFD), we can now develop improved design tools for 
rocket propulsion. However, there are three requirements that the CFD tool must meet before it can be routinely employed 
in the design process. First, there must be an increased fidelity of the solution over the historical tools in terms of numerical 
accuracy and physical modeling. Second, a rigorous validation process is required before the designer can confidently 
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employ CFD as a tool for quantitative decision-making. Finally, a significant quantity of parametric solutions must be 
obtainable during the design cycle time frame. 

The application of CFD to injector design has lagged behind other areas such as turbomachinery because the physical 
problems are more complicated for multiphase, turbulent reacting flows. Although many challenges remain, it is our view 
that with concerted community efforts, significant gain in applying computational tools for the design of injector and other 
combustion devices can be realized in the near future. For example, new models that efficiently account for some of the 
complex processes^ and, thus, increase the solution fidelity, have recently become available. The three-dimensional 
geometry of multi-clement injectors and the complex physical processes inherent in the flows that issue from them create 
major problems in the other two areas in terms of using CFD as a design tool. The harsh high pressure and temperature 
environments typical of injector flows create significant difficulty in obtaining experimental data of satisfactory quality to 
validate and guide further development of computational models. Additionally, solving the equations, for multiphase 
reacting flows, with high resolution typically requires long computational time. However, the continued increase in 
computer speed and progress in parallel processing has begun to soften this turnaround problem. 

It has long been known that small changes in injector geometry can have significant impact on performance^, as well 
as on environmental variables such as combustion chamber wall and injector face temperatures and heat fluxes. CFD can be 
fruitfully used to account for several independent geometrical variables over appropriate ranges to assess the impact on 
design goals. To effectively accomplish this multi-objective exercise, an optimization technique is desirable to assist and 
guide the designer in objective management of the complex and competing trends embodied in the CFD-generated data. 


2. SCOPE & APPROACH 

In the present effort, a CFD-based injector design/optimization methodology for single element injectors is proposed 
and applied to an injector targeted for the future reusable launch vehicle technologies. While the ultimate goal is to analyze 
multi-element injectors, much of the detailed work in injector design can be done, or at least initiated, at the single element 
level. In the present effort, the focus is on the single element methodology. CFD solutions based on the Navier-Stokes 
equations, finite rate chemistry, and the k-e turbulence closure are generated with design of experiment (DOE) 
techniques^’*®, and the response surface (RS) method’ is employed as the optimization tool. The design considerations are 
guided by four design objectives, namely, the maximum temperature on the oxidizer post tip, the maximum temperature on 
the injector face, the adiabatic wall temperature, and the length of the combustion zone. Together, they address the concerns 
arising from injector performance and life. Four design variables are selected, namely, H2 flow angle, H2 and O2 flow areas 
with fixed mass flow rates in fuel and oxidizer, and O2 post tip thickness. As a proof of the concept, the present CFD-RS 
optimization strategy is demonstrated in the context of a coaxial element using gaseous oxygen (GO2) and gaseous 
hydrogen (GH2) propellants. The response surface generation is evaluated for two grids of different resolutions and the 
relationship between the optimization process and CFD solution process is explored. Finally, the response surfaces are used 
in three separate optimization examples demonstrating the utility and flexibility of the methodology. 

To facilitate the development of the present methodology, a baseline element design is needed as a starting point. This 
baseline concept is generated by an empirical design methodology based on a specific set of propellant flow rates, mixture 
ratio and chamber pressure. Dependent variables to characterize performance and environments can then be estimated from 
this design. The selected design variables are then varied based on this baseline design and the design space populated with 
the aid of DOE technique. The CFD cases are run and post processed to extract the required dependent variable data. This 
data is then used to generate a response surface for each dependent variable in terms of the independent variables, and an 
optimization technique is applied to identify suitable designs. Details on these ideas/processes will be provided after 
describing the injector model. 


Injector Model 

In liquid rocket propulsion, there are two basic types of injector elements. The first type is an impinging element 
(Figure la) where mixing occurs by direct impingement of the propellant streams at an acute angle. The impingement 
enhances mixing by head-on interaction between the oxidizer and fuel®. The second type of injector consists of non- 
impinging elements where the propellant streams flow in parallel, typically in coaxial fashion (Figure lb). Here, mixing is 
enhanced through shear-mixing processes**. 

From a design standpoint, both element types have some appealing characteristics. For instance, if the impinging 
element has an F-O-F arrangement, the mixing occurs rapidly, which can yield high performance. However, since the 
combustion zone is close to the injector face, the potential for high levels of injector face heating must be addressed. If the 
non-impinging element is assumed to be a shear coaxial element, mixing across the shear layer is relatively slow®, requiring 
longer chambers to allow complete combustion. Since the combustion zone is spread over a longer axial distance, the 
injector face is generally exposed to less severe thermal environments. Important parameters for the impinging element 
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(assuming fixed mass flow rates and constant propellant inlet conditions) include relative orifice size (or, relative stream 
momentum ratio), impingement angle and orifice spacing. Important parameters for the shear coaxial element (assuming 
fixed mass flow rates and constant propellant inlet conditions) include the area ratio of the two concentric tubes (or velocity 
ratio) and the shear area between the two propellant streams (the oxidizer post tip thickness)**. 

Of course, one can combine the above-mentioned characteristics of injector types to develop hybrid concepts. For 
example, it has been noted* that modifications for performance improvement in the shear coaxial element can be realized by 
directing the fuel toward the oxidizer stream rather than parallel to it and thinning the oxidizer post wall. Both 
modifications allow the shear coaxial element to take on some of the aspects of the F-O-F impinging element. These 
notions lead to the hybrid element shown in Figure 2 that has been developed by The Boeing Company (U. S. Patent 
6253539). 

Four independent design variables are chosen for the element shown in Figure 3a, namely, the angle at which the H 2 is 
directed toward the oxidizer, Variable- 1 (a), the change in H 2 flow area from the baseline, Variable-2 (AHA), the change in 
O 2 flow area from the baseline, Variable-3 (AOA), and the oxidizer post tip thickness, Variable-4 (OPTT). The fuel and 
oxidizer flow rates are held unchanged. The ranges of the independent variables are shown in Table 1. The dependent 
variables chosen are the maximum face temperature. Objective- 1 (TFmax)» the wall temperature at a distance three inches 
from the injector face, Objective-2 (TW4), the maximum oxidizer post tip temperature, Objective-3 (TTmax) and centerline 
axial location where the combustion is 99% complete, Objective-4 (hereafter referred to as the combustion length, X,,) 
(Figure 3b). The combustion length, Objective-4 (Xcc), was chosen as a measure of performance. Shorter combustion 
lengths indicate better performing designs. The three temperatures (calculated as adiabatic wall temperatures in this study) 
were chosen as indicators of local environments. Lower temperatures would indicate a design that had longer life due to 
decreased thermal strain on the part. 

Numerical Procedure 

A pressure-based, finite difference, Navier-Stokes solver, FDNS500-CVS*^'*'*, is used in this study. The Navier-Stokes 
equations, the two-equation turbulence model, and kinetic equations are solved. Convection terms are discretized using 
either a second order upwind, third order upwind or a central difference scheme, with adaptively added second order and 
fourth order dissipation terms. For the viscous and source terms, second order central differencing scheme is used. First 
order upwind scheme is used for scalar quantities like turbulence kinetic energy and species mass fractions, so as to ensure 
positive values. Steady state is assumed and an implicit Euler time marching scheme is used for computational efficiency. 
The chemical species continuity equations represent the H 2 - O 2 chemistry. It is represented with the aid of a 7-species and 
9-reaction set*^ *^. The simulation domain and the boundary conditions used in all the CFD cases are shown in Figure 4. 
Because of the very large aspect ratio, both the injector and chamber have been shortened (at the cross hatched areas) for 
clarity. Both fiiel and oxidizer flow in through the west boundary where the mass flow rate is fixed for both streams. The 
nozzle exit, at the east boundary, is modeled by an outlet boundary condition. The south boundary is modeled with the 
symmetry condition. All walls (both sides of the oxygen post, the outside of the fuel annulus, the outside chamber wall, and 
the faceplate) are modeled with the no slip adiabatic wall boundary condition. 

Desien of Experiments 

The predictive capability of a response surface (RS) is largely dependent on the selection of the design space as it 
dictates the distribution of the available information. There are different types of design of experiments (DOE) techniques. 
Among various alternatives, the orthogonal array (OA)*^, which is an efficient approach, is used here. An OA is a fractional 
factorial matrix that assures a balanced comparison of levels of any factor or interaction of factors*^. Consider A, a matrix 
with elements of a{ where j shows the row (/ = 1,2... and i shows the column (i = 1,2. ../i^) that a/ belongs to, 
supposing that each aj eQ = {0,1...^-! }. A is called an orthogonal array of strength t <ric if in each n^row-by-r-column 

sub-matrix of all possible distinct rows occur A times. Such an array is denoted by OA(nryfic,q,t) by Owen*^. Based on the 
DOE theory, orthogonal arrays can significantly reduce the number of experimental configurations. Information directly 
relevant to the use of OA can be found in the works of Shyy et al.*^ and Papila*^. 

In the present study, using OA 54 designs are generated for fitting and testing the response surface. To test the RS, 14 
of the 54 designs are selected using cross-validation techniques. Cross-validation*^ is an established technique for 
estimating the prediction accuracy. This method is usually performed using either a number of random test/train partitions 
of the data, or ^-fold cross-validation**. During the CFD computations, 2 of the 40 designs (training set) were found to be 
unacceptable because they exhibited unsteady behaviors while the numerical algorithm was based on the steady state 
model. Hence the final information included 38 designs for training the RS and 14 to test their predictive capabilities. All 
the design variables are normalized between 0 and 1 based on their maximum and minimum values. All the responses 
obtained from the CFD solutions of the 52 valid designs are scaled to 0(1) and RS generated. Once the RS are generated, 
the responses are renormalized between 0 and 1 based on the maximum and minimum of the generated RS, which will be 
referred to as the normalized responses in the text. Henceforth these normalized values will be used in this paper. The 
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training and testing designs are shown in Table A1 in the appendix. It should be noted that when the Variable-3 (AOA) is 1 
or 0, the O 2 flow area is reduced by 0% or 40%, respectively, as compared to the baseline area. 

Response Surface Methodology 

The response surface methodology (RSM) is adopted in the present study. This technique is effective in representing 
the global characteristics of the design space and it filters noise associated with design data. However, depending on the 
order of polynomial employed and the shape of the actual response surface, the RSM can introduce substantial errors due to 
its generally low order representation. However, the number of design points needed to construct a response surface 
increases rapidly with the number of design variables. For example, for a second order response surface, the number of 
coefficients associated with the polynomial scales with the number of design variables, N, as (N-\-l)(N-\-2)/(2!), while a third 
order polynomial scales as (N-\-lXN+2)(N-\-3)/(3!), 

In this study the response surfaces are fit with n*** order polynomials using standard least-squares regression. Statistical 
analysis software, is used for the generation of the polynomials. In the standard least-squares regression method, 

based on t-statistics^ terms with least influence in the prediction are discarded and the prediction improved, which gives the 
reduced model. The quality of fit between different surfaces can be evaluated by comparing the adjusted RMS error^, Oa, 
defined as: 



( 1 ) 


where is the error at the point, n is the number of data points and p is the number of coefficients. The measure of error 
given by is normalized to account for the degrees of freedom in the model. This adjusted RMS error thus accounts for 
the nominal effect of higher order terms providing a better overall comparison among the different surface fits. Another 
measure of error, the adjusted coefficient of multiple determination^, /?/, measures the proportion of variation in the 
response around the mean that can be attributed to terms in the model rather than to random error. It is given by 
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where SSyy is the sum of squares of error about the mean of the observations. For a good fit, Ra should be closer to 1 . The 
RMS error, <r, for the test data is given by: 



(3) 


In this equation Ci is the error at the test point and m is the number of test points. For a good model the difference 
between Ga and o*should not be large. 

In certain cases, the number of CFD computations may not be adequate enough to spare enough data for testing the RS. 
Hence an alternate method to estimate the performance of the RS is to compute the PRESS statistic. The method proposed 
by Allen^®’ compute a residual sum of squares. This residual is obtained by fitting a RS over the design space after 
dropping one design point from the training set and then comparing the RS predicted value for that point with the expected 
value. The PRESS rms error is given by 
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where is the expected value, p, is the value predicted by the RS for the i* point which is excluded while generating the 
RS and n is the number of design points. If this value is close to then the model performs well. 

Construction of Multi-Obiective Optimization 

The optimization problem at hand can be formulated as 
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X = (Xi, X2, . . x„)'^| 

subject to lb<\\< ub. 


where lb and ub are the lower and upper bound, respectively, of the design variable Xj. In the current study, different 
combinations of the objective functions are used to investigate the issues related to multi -objective optimization. One 
method of optimizing multiple responses simultaneously is to build, from the individual responses, a composite response 
known as the desirability function , The method allows for a designer’s own priorities, on the response values, to be built 
into the optimization procedure. The first step in the method is to normalized the responses by develop a desirability 
function, di for each response, i. In the case where a response is to be minimized the desirability takes on the form: 


= 


f-a Y 

B-A 


,B<F<A 


(5) 


where F is the prediction of the RS, F is the target value and A is the highest acceptable value such that = I for any F <B 
and di = 0 for F > A. Hence to minimize a response, this function has to be maximized. Choices for A and B are made 
according to the designer’s priorities (i.e. maximum and minimum, respectively, of the individual RS, in the present study). 
The sensitivity of the parameter, 5 , illustrated in Figure 5 can be instructive’. Desirability functions with s« I imply that a 
product need not be close to the response target value, R, to be quite acceptable. But 5 = 8, say, implies that the product is 
nearly unacceptable unless the response is close to B. 

Value of s is set based on which response has higher priority. A single composite response is developed which is the 
geometric mean of the desirability functions of the individual responses. The composite response is defined as: 

D = {d,d2-d3...jyi ( 6 ) 


where m is the number of objectives, which is then maximized. This approach was adopted in our previous efforts*^ and has 
been found to be instructive. 

The optimization problem is then solved using Solver, an optimization tool available as part of Microsoft Excel 
package^^. This tool uses the Generalized Reduced Gradient {GRG2) nonlinear optimization code developed by Lasdon et 
al.^^ The optimum solutions are then compared with CFD computations. Since the responses are scaled to 0(1), the error 
measures have to be accordingly scaled to estimate the accuracy of the obtained solutions. In terms of the actual values the 
error for a response is defined as 


\yCFD 

error = - 
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where ycFo is the solution obtained from the CFD computations and y^s is the prediction of the RS. Using simple 
mathematics, not shown here, the error in the scaled variables can be written as 


error = 


ycFD yRs 

yCFD ^ 


( 8 ) 


where the bar represents the scaled values, and K is defined as 




y 


max 


ymin 


(9) 


where y^nm and y^ are the actual minimum and maximum values, respectively, based on the available set of training and 
testing data for that response. 

3. RESULTS AND DISCUSSION 

Grid Sensitivity Investigation 

Initially the 54 cases identified by DOE were run on an axisymmetric geometry with 336x81 nodes. Only 33 out of the 
40 cases gave valid results. The rest of the results showed unsteady features, which cannot be properly captured by the 
steady state model employed. The RSs generated with these data for Objective-3 (TTmax) and Objective-4 (Xcc) had Ra^ 
values of 0.961 and 0.976, respectively, suggesting a moderate fit. On checking the grid distribution in the combustion zone 
for cases that gave us steady state solutions, it was determined that the grid resolution was insufficient. After a series of 
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tests involving addition of grid points in the axial direction in the combustion zone, a 430x81 grid was deemed appropriate 
and used for the second run of the optimization study. To highlight the grid refinement, a comparison of the grid 
distributions is shown in Figure 6. The final grid was the product of tripling the axial node density in the combustion zone. 
The thick lines show the initial grid density, while the thin lines show final grid density. Note that, for clarity, only every 
sixth j-line is shown. The RSs for Objective-3 (TTn^jx) and Objective-4 (Xcc) had Ra^ values of 0.989 and 0.995, 
respectively, representing a considerable improvement over the RS performance based on the first, coarser grid. It is to be 
noted that the Rg^ values for Objective- 1 (TF^ax) and Objective-2 (TW4) are 0.999 for both the initial and final grid. This 
time, only 2 out of the 40 designs failed to provide valid results. This experience indicates that in addition to facilitating 
design optimization, the RSM can also help address the adequacy of the CFD solution accuracy. It offers insight into 
potential problems, based on the statistical regressions, based on which we can refine the computations and improve the 
fidelity of the individual and collective databases. Of course, this approach assumes that the given order of the RS assigned 
is adequate for representing the design space. While this approach by no means guarantees universally satisfactory 
outcomes, it suggests clear directions to assess the accuracy improvement aspect of data generation, which can be critical. 

Interplay of Design Variables and Impact on Design 

CFD computations are done for the 38 cases listed in Table A1 and the 4 objectives are estimated. On observing the 
solutions, it is noted that the flow field is dominated by a large recirculation zone (Figure 7), This recirculation zone 
occupies the bulk of the downstream domain so as to accommodate the sudden expansion of the geometry. There is also a 
smaller circulation zone noticed adjacent to the injector. 

A comparison between two particular Cases 12 and 44, helps demonstrate the usefulness of the proposed methodology. 
Table 2 shows the independent variable values and it is clear that values of the design variables vary prominently for these 
two cases. Normalized results for both cases, including the combustion length are also shown in Table 2. The temperature 
contours for the two cases are shown in Figures 8 and 9. Results from Case 44 are shown on top and Case 12 on bottom in 
both figures. Figure 8 shows that the temperatures on the combustor wall and injector face for Case 44 are lower than those 
for Case 12 indicating desirable attributes from the life consideration. On the other hand, in the near injector region shown 
in Figure 9, the post tip temperature for Case 44 is higher than that for Case 12, suggesting that there is no effective 
improvement of life for Case 44. Case 12 has a short combustion length, suggesting favorable performance, but has very 
high injector face and wall temperatures, reducing an injector’s life expectancy. In Case 44, where the injector face and 
wall temperatures are low, a longer combustion length affects the performance. These conflicting observations in the 
injector details can have substantial effects on the design outcome. These observations suggest that Objective- 1 (TF^ax) and 
Objective-2 (TW4) that reflect the life of the element follow a similar trend. This can be attributed to the larger, outer 
recirculation zone (Figure 7), which strips off hot combustion gases from the flame region, recirculates them back along the 
wall and injector face, and hence regulates Objective- 1 (TFnux) and Objective-2 (TW4). However, the third metric 
representing life follows a reverse trend suggesting that it is tough to obtain a design, which would be able to improve the 
life of the element. Also, it is apparent that neither of these designs can accommodate both performance and life at the same 
time. An optimization technique will be needed to guide the designer, through the competing trends exhibited by the 
dependent variables, to reach a logical compromise. More discussion will be made when inspecting design optimization 
tasks based on individual objectives. 

Optimization Process 

In the following, the 38 acceptable data are used to generate the RS for the 4 objective functions. Both full and reduced 
quadratic polynomials are generated. Table 3 identifies Ga and a for the different objective functions and different 
polynomials generated. The full quadratic model is consistent in performance in terms of both Ga and G. The reduced 
quadratic models have either a poorer value of Gox offer only marginal improvement over the full response surface model. 
These observations suggest that there can be modeling deficiency and/or noise in the data. 

Now, considering the full quadratic models of the various objectives, the fit for Objective- 1 (TFmax)» Objective-2 (TW4) 
and Objective-4 (Xcc) is satisfactory, suggesting that there is no need for further improvement. However, while comparing 
the full quadratic RS predictions to the CFD results a spread is noticed in the plot of Objective-3 (TTmax) (Figure 10). 

Based on this observation, a cubic model is generated for Objective-3 (TTmax)- A cubic model in 4-design variable 
needs 35 design points. To obtain a good fit the number of design points should be considerably larger than the required 
number. Since there are only 38 design points available from the training set, the testing set is also included in the training 
set and all 52 design points are used to generate a reduced cubic model. The PRESSnm is used to estimate the performance 
of the generated RMS along with the Table 4 compares the performance of the quadratic and the reduced cubic models. 
The reduced cubic model is seen to perform better than the quadratic model. Hence this cubic model is used for Objective-3 
(TTmax) in the optimization phase. The RSs are listed in the Appendix. 

Using these RSs, an optimization study is conducted to study the influence of objectives controlling life to the ones 
controlling performance and also to confirm the observations made before based on the CFD solutions. First each of the 
objectives is minimized individually (Table 5). Following this, the influence of thermal environments on the performance is 
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estimated (Table 6 ). Next, composite optimization is performed with different weights emphasizing either the performance 
or life of the component (Table 7). 

{a) Single-Objective Optimization 

Figure 11 shows the responses and the optimum designs listed in Table S. Since the chemical heat release is spread all 
over the flow field, with the aid of the large recirculation zone (Figure 7), Objective- 1 (TF^ax) and Objective-2 (TW4) 
correlate with each other directly (Opt-Cases 1 and 2). From this study it is clearly seen that when one of these two 
objectives is minimized, the other objective is low as well. However, the third metric of life, Objective-3 (TTmax) is high for 
these cases (Figure 11a). Due to conflicting trends of the temperatures, it seems difficult to obtain an optimum design based 
on all three life objectives. It is also noticed that when Objective- 1 (TFo^x) and Objective-2 (TW4) are low, combustion 
length is long (Opt-Cases 1 and 2). Figure 11b shows that both Objective-1 (TFnu«) and Objective-2 (TW4) are minimized 
for low H 2 flow angle and thick O 2 post tip. This results in the propellants flowing in parallel just as in a shear coaxial 
injector, resulting in slower mixing across the shear layer and less severe thermal environment. Table 5 and Figure lib 
indicate that Variable-3 (AO A) exhibits different trends between Objective- 1 (TF^ax) and Objective-2 (TW4), and between 
Objective-3 (TTmax) and Objective-4 (Xcc), while the other three design variables behave consistently between the two 
groups of the objectives. As shown in Figure 12a Objective-1 (TFmax) is less sensitive to the oxidizer flow area (AOA) than 
Objective-2 (TW4). Overall, a small oxidizer flow area (Variable-3 (AOA) = 1) minimizes Objective-2 (TW4). On the other 
hand, large fuel flow area, minimizes both Objective-1 (TFmax) and Objective-2 (TW4) (Figure 11b). Figure 12b shows that 
Variable-3 (AOA) has opposing influence on Objective-3 (TTmax) and Objective-4 (Xcc). Objective-3 (TTmax) is minimized 
for large oxidizer flow area whereas Objective-4 (Xcc) is minimized for small oxidizer flow area. Both Objective-3 (TTmax) 
and Objective-4 (Xcc) are minimized for large H 2 flow angles, small fuel flow area and thin O 2 post tip. 

(b) Multi-Objective Optimization 

To accommodate life and performance considerations, a multi-objective optimization study is carried out. Starting with 
Objective-4 (Xcc) (performance), the objectives influencing thermal environments are added one at a time. The obtained 
optimum designs are presented in Table 6 and Figure 13. When Objective-3 (TTmax) is minimized along with Objective-4 
(Xcc) (Opt-Case 5), the optimum design seems to favor the later (see Figure 1 lb and Table 5). More insight can be gained 
by visualizing Figure 12b. The individual objectives drive the design towards the opposite ends of the range of Variable-3 
(AOA). But the composite desirability function has marginal variation to Variable-3 (AOA). Hence, any particular value of 
Variable-3 (AOA) will not cause a large impact, and the optimum design is affected mostly by other design variables. 

When Objective-1 (TFmax) is minimized along with Objective-4 (Xcc) and Objective-3 (TTmax) (Opt-Case 6 ), a 
moderate drop in the performance is noticed. The compromise design has a high H 2 flow angle, and thin O 2 post tip, which 
favors Objective-3 (TTmax) and Objective-4 (Xcc) (see Figure 11b and Table 5). However, the H 2 flow area, Variable-2 
(AHA), is large, which favors Objective-1 (TFmax) (see Figure 11b and Table 5). Addition of Objective-2 (TW4) (Opt-Case 
7) in the minimization process does not alter the optimum design. Figure 14a, shows clearly the reason why the Variable-2 
(AHA) reaches the maximum in the present multi-objective optimization (Opt-Cases 6 and 7). Variations of Objective ! 
(TFmax) and Objective-2 (TW 4 ) to Variable-2 (AHA) are high as compared to the other two objectives. This drives the 
composite desirability function in favor of larger H 2 flow area and hence lowers temperatures at the injector face and the 
adiabatic wall without severely degrading the performance. Figure 14b shows that temperature metrics have marginal 
variation to the oxidizer flow area whereas the performance (Objective-4 (Xcc)) is minimized for low values of the variable. 
Hence the composite optimization procedure is driven in favor of (Objective-4 (Xcc) (see Figure 1 lb and Table 5). 

To further probe the interplay between performance and life, a composite optimization study is conducted with varying 
weights on the objective functions. Figure 15a shows the values of the objectives for the optimum designs listed in Table 7. 
Opt-Cases 8 and 10 return almost the same optimum designs (Figure 15b). The optimum design also agrees with Opt-Case 
7. This shows that the optimum design is reasonably robust with respect to moderate variations in weighting parameters. In 
Opt-Case 9 where the temperatures are weighted considerably heavier than Objective-4 (Xcc), the optimum design 
compromises between the designs obtained for Opt-Cases 1, 2 and 3. The H 2 flow angle is small making the propellant flow 
parallel to each other. At the same time, the O 2 post tip is thin which lowers Objective-3 (TTn^x) but increases the shear 
area between the propellant streams, thereby increasing the mixing rate. Hence, Objective-3 (TTnuw) is not lowered to the 
same extent as the other temperatures. The H 2 flow area favors Objective-1 (TFmax) and Objective-2 (TW 4 ). Due to smaller 
weight assigned to Objective-4 (Xcc), the oxidizer flow area is not small as in Opt-Case 7. Figure 16 shows that as the 
oxidizer flow area increases, performance (Objective-4 (Xcc) decreases. The variation of the temperatures is moderate to the 
oxidizer flow area. The composite desirability function which is formed with large weights on the temperatures has a 
marginal variation to the flow area as well. When Objective-4 (Xcc) is weighted considerably heavier (Opt-Case 1 1), the 
temperatures (Objective- 1 (TFmax) and Objective-2 (TW 4 )) turn out to be high. This affects the life of the injector 
noticeably. The optimum design has a moderate H 2 flow angle and thin O 2 post tip favoring performance. The results show 
that in cases where Objective-4 (Xcc) is emphasized (Opt-Cases 8 , 10 and 11), the optimum O 2 flow area is always low 
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hence the velocity is high (Figure 15b). This results in higher mixing of the oxidizer with the fuel, when the H2 stream is 
aimed inward towards the O2 stream. 

CFD solutions have validated the optimum designs obtained from RS for all cases. The error, as defined by Eq. 8, is 
marginal. The error, using Eq. 8, is to the order of about 3-4% in the case of Objective-3 (TTmax) and to the order of 1% for 
the other objectives. The discrepancy in Objective-3 (TTmax) is most probably due to the steady state assumption for an 
unsteady problem. Opt-Case 2 shows some inconsistency in the values of Objective-4 (Xcc), which is due to the unsteady 
features noticed in the solution. 

4. SUMMARY AND CONCLUSION 

In this study a demonstration of a CFD-based design optimization methodology for a single element rocket injector is 
presented. The design variables controlling the thermal environment and performance of the element are identified and the 
design points selected using DOE. CFD computations are carried out for each design and appropriate orders of polynomial- 
based response surfaces are generated for the estimated objectives. Based on the performance of the RSs a grid refinement 
study is carried out and a new efficient grid is generated. Following the results obtained from the CFD computations, an 
optimization study is carried out and correlations between different objectives are identified. Relationship between the dual 
goals of performance and life is studied elaborately. Shortcomings due to the steady state approximations are noticed, 
suggesting the requirement of developing a robust unsteady CFD model. Based on the results obtained we have reached the 
following conclusions. 

• A large circulation in the combustion chamber distributes the heat flux over the flow field which influences the 
temperature on the injector face and the adiabatic wall of the combustor. 

• Objective-1 (TFmax) and Objective-2 (TW4) are minimized for small H2 flow angle and thick O2 post tip. The thick O2 
post tip increases die area for heat transfer and hence, Objective-3 (TTmax) is high. Additionally, Objective-3 (TTmax) is 
minimized for a design which maximizes Objective-1 (TFmax) and Objective-2 (TW4). This conflicting trend makes the 
optimization based on life a challenge. When all the life metrics are given importance, a design similar to a shear 
coaxial injector is obtained with thin O2 post tip, 

• Optimization of the injector design based on Objective-1 (TFmax) or Objective-2 (TW4), results in a design similar to a 
shear coaxial injector, whereas the one based on Objective-4 (Xcc) (performance metric) results in a design similar to 
an impinging injector. 

• Optimization study provides insight into the interplay of the design variables and objectives. For example, when 
individually optimized, Objective-3 (TTmax) and Objective-4 (Xcc) responds to Variable-3 (AOA) differently. But the 
sensitivity of the desirability function, combining the two objectives, to Variable-3 (AOA) is marginal, allowing the 
design to assume flexible value of Variable-3 (AOA) without causing large impact. Sensitivity of Objective- 1 (TFmax) 
and Objective-2 (TW 4 ) to Variable-2 (AHA) is high as compared to the other two objectives, which drives the 
composite desirability function based on all four objectives in favor of larger H2 flow area without severely degrading 
the performance. When the optimization process is biased towards either the performance or the life metrics, it 
deteriorates the other. For example, when the life metrics are weighed heavily a compromise design is obtained which 
is similar to a shear coaxial injector with a thin O2 post tip that has a long combustion length. 

• It seems that the best compromise is obtained by weighing the life and performance metrics equally. It results in a 
design which is similar to an impinging injector with thin O2 post tip. 

• Finally, RS can aid in judging the quality of the CFD solution and help improve modeling issues, for e.g. improving 
grid quality. 
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Variable 

Minimum 

Maximum 

Actual 

Normalized 

Actual 

Normalized 

a 


0 


1 

AHA 

fBaseHne 

0 


1 

AOA 

Baseline-40% 

0 

Baselihe 

1 

OPTT 

O.Olin 

0 

/ 0.02 ip 

1 


( J Table 1: Range of design variables .V 1 , 

4o re a/.o 'J ^ d 


Case 

a 

AHA 

AOA 

OPTT 

TFma, 

TW 4 

TTmax 

Xcc 

12 

1 

0 

0 

0 

0.998 

0.927 

0.128 

-0.004 

44 

0 

0.5 

0.5 

1 

0.285 

0.395 

0.923 

0.567 


Table 2. Independent and dependent variable (Objectives) for Cases 12 and 44 from CFD computations 


{Normalized values shown}. 
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Full quadratic RS 

Reduced quadratic model 

max 

Number of observations 

38 

38 

a. 

0.00566 

0.00546 

a (14 points) 

0.00460 

0.00463 

Mean 

0.495 

TW4 

Number of observations 

38 

38 

Ca 

0.00803 

0.00795 

a (14 points) 

0.00669 

0.00799 

Mean 

0.514 

T^max 

Number of observations 

38 

38 

Oa 

0.0413 

0.0401 

a (14 points) 

0.0396 

0.0382 

Mean 

0.560 

Xcc 

Number of observations 

38 

38 

Oa 

0.0205 

0.0197 

a (14 points) 

0.0178 

0.0186 

Mean 

0.497 


Table 3: Performance of full and reduced quadratic RSs for the four objective functions {Non-normalized values used } . 




Full quadratic RS 

Reduced cubic RS 

T^max 

Number of observations 

38 

52 


Oa 

0.0413 

0.0303 


PRESS, „„ 

0.0521 

0.0388 


Mean 

0.560 

0.591 

Table 4: Performance of RSs for the TTma^- Reduced cubic RS has 21 coefficients 

[Non-normalized values used). 


Opt-Case 

a 

AHA 

AOA 

OPTT 

TF„a, 

TW4 

TTmax 

Xcc 

1 

0 

1 

0.592 

1 

0.0(1) 

0.0725 (0) 

0.914 (0) 

0.769 (0) 

CFD 





-0.00207 



0.758 

error (%) 





0.11 



0.25 

2 

0 

1 

0 

1 

0.0309 (0) 

0.0(1) 

1.0(0) 

0.440 (0) 

CFD 





0.091 

0.0461 

0.911 

0.568 

error {%) 






0.57 

2.84 

2.93 

3 

1 

0 

1 

0 

0.944 (0) 

0.976 (0) 

0.0(1) 

0.153 (0) 

CFD 





0.943 

0.969 

0.103 

0.158 

error (%) 





0.01 


4.46 

0.12 

4 


0 

0 

0 

0.987 (0) 



0.0(1) 

CFD 





0.981 

0.919 

0.119 

-0.004 

error (%) 





0.08 


2.69 

0.12 

Table 5: Minimizing individual 0 

)jectives. Value in parent 

hesis (1) indicates which objective function is minimized 


{Normalized values shown}. 
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Opt-Case 

a 

AHA 

AOA 



TW4 


Xc. 

4 

0.917 

0 

0 

0 

0.987 (0) 

0.926 (0) 

0.182 (0) 

0.0(1) 

CFD 





0.981 

0.919 

0.119 

-0.004 

error (%) 





0.08 

0.08 

2.69 

0.12 

5 

1 

0 

0 

0 



0.151 (1) 

0.00097 (1) 

CFD 





0.998 

0.927 

0.128 


error (%) 





0.03 



0.14 

6 

1 

1 

0 

0 

0.376(1) 



WililKfellWi 

CFD 





0.369 




error (%) 





0.10 


3.93 

0.38 

7 

1 

1 

0 

0 

0.376(1) 

0.224(1) 


WWrEfiMI 

CFD 





0.369 

0.214 

0.252 

0.264 

error (%) 





0.10 

0.12 

3.93 

0.38 


Table 6: Study of effect of life on performance. Value in parenthesis indicates the weighting on the objective functions 

{Normalized values shown}. 


Opt-Case 

a 

AHA 

AOA 

OPTT 

■EillHlIi 

WMSSBKH 

TT„„, 

Xcc 

8 

1 

1 

0.192 

0 

0.346(1) 

0.210(1) 

0.163(1) 

0.334 (0.5) 

CFD 









error (%) 









9 

0 

1 

0.497 

0 

0.0451 (5) 

0.0823 (5) 

0.471 (5) 


CFD 









error (%) 









10 

1 

1 

0 

0 

0.376 (0.5) 

0.224 (0.5) 

0.155 (0.5) 

0.279(1) 

CFD 





0.369 

0.214 

0.252 

0.264 

error (%) 






0.12 

3.93 

0.38 

11 

0.612 

0 

0 

0 

0.919(0.1) 

0.898 (0.1) 

0.281 (0.1) 


CFD 









error (%) 







1 



Table 7: Study of influence of life and performance objectives on each other. Value in parenthesis indicates the weighting 

on the objective functions (Normalized values shown). 
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Figure 4: Simulation domain and boundary conditions 
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Figure 6; Comparing the unrefined (336x81) {thicker lines) and refined (430x81) (thinner lines) grids. 
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Response _ Objective 


Optmuzation Study 


D TFmax 
TV 

o TTma 
Xc( 


X OPTT 


Opt-Case 


2 3 

Opt-Case 


(a) (b) 

Figure 11 : Minimization of different objectives. (Case 1: TFroa*, Case 2: TW 4 , Case 3: TTnax. Case 4; X„) 
(a) Objectives, (b) Design Variables {Normalized values shown}. 


Response of TFmax and TW4 


Response Functions 



TFmax 


g, 0.5 



Figure 12: Variation of (a) TFra« and TW 4 for o=0, AHA=1 and OPTT=l, (b) TTnax and X«c for ct=l, AHA=0 and 
OPTT=0, with respect to AOA (Normalized values shown and D is the desirability function}. 












(a) (b) 

Figure 13: Composite minimization of objectives with different weightings.(Case 4: (0,0,0,1), Case 5: (0, 0 , 1 , 1 ), Case 6 : 
(1,0,1, 1), Case 7: (1,1,1,!),) {The values in parenthesis indicates weights for (TF^ax. TW 4 , TTmax. Xcc)), 

(a) Objectives, (b) Design Variables (Normalized values shown). 



Response Functions 



AOA 


(b) 

Figure 14: Variation of (a) TFnu,, TW4, TTroa and Xcc, with respect to AHA for o=l, AOA=0 and OPTT=0, (b) TF^ax. TW 4 . 
TTmax and Xcc, with respect to AOA for ot=l, AHA=1 and OPTT=0 (Normalized values shown and D is the desirability 

function). 
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Opt-Case 


(a) (b) 

Figure 15: Composite minimization of objectives with different weightings.(Case 8: (1,1, 1,0.5), Case 9: (5,5,5,0.1), Case 
10: (0.5,0.5,0.5,1), Case 11: (0.1,0. 1,0. 1,5)) {The values in parenthesis indicates weights for (TF^u, TW4, TTroj„, Xcc)}, 
(a) Objectives, (b) Design Variables (Normalized values shown). 



Figure 16: Variation of TFa^j,, TW4, TTnux and X«c, with respect to AOA for 0=0, AHA=1 and OPTT=0 (Normalized 

values shown and D is the desirability function). 




APPENDIX 



H 2 flow angle (cd 

Change in H 2 flow 
Area (AHA) 

Change in O 2 flow 
Area (AOA) 

O 2 post tip thickness 
(OPTT) 

Case 1: 

0 

0 

1 

0 

Case 2: 

0 

0 

1 

0.5 

Case 3: 

0 

0 

1 

1 

Case 4: 

0.5 

0.5 

0 

1 

Case 5: 

1 

1 

0.5 

0 

Case 6: 

1 

1 

0.5 

0.5 

Case?: 

1 

1 

0.5 

1 

Case 8: 

0 

0.5 

0.5 

0 

Case 9: 

0.5 

1 

1 

0 

Case 10: 

0.5 

1 

1 

0.5 

Case 11: 

0.5 

1 

1 

1 

Case 12: 

1 

0 

0 


Case 13: 

1 

0 

0 

0.5 

Case 14: 

1 

0 

0 

1 

Case 15: 

0 

1 

0 

0 

Case 16: 

0 

1 

0 

0.5 

Case 18: 

0.5 

0 

0.5 

0 

Case 19: 

0.5 

0 

0.5 

1 

Case 20: 

1 

0.5 

1 

0 

Case 21: 

1 


1 

0.5 

Case 22: 

1 

0.5 

1 

1 

Case 23: 

0 

1 

0.5 

0 

Case 24: 

0 

1 

0.5 

1 

Case 25: 

0.5 

0 

1 


Case 26: 

0.5 

0 

1 

1 

Case 27: 

1 

0.5 

0 

0 

Case 28: 

1 

0.5 

0 

1 

Case 29: 

0 

0 

0 

0 

Case 30: 

0 

0 

0 

0.5 

Case 32: 

0.5 

0.5 

0.5 


Case 33: 

1 

1 

1 

0 

Case 34: 

1 

1 

1 

1 

Case 35: 

0 

0.5 

1 


Case 36: 

0 

0.5 

1 

1 

Case 37: 

0.5 

1 

0 

0 

Case 38: 

0.5 

1 

0 

1 

Case 39: 

1 

0 

0.5 

0 

Case 40: 

1 

0 

0.5 

1 


Table Ala: Training data (missing Case numbers indicate the unacceptable Cases). 
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a 

AHA 

AOA 

OPTT 

Case 41: 

0.5 

0.5 

0 

0 

Case 42: 

0.5 

0.5 

0 

0.5 

Case 43: 

0 

0.5 

0.5 

0.5 

Case 44: 

0 

0.5 

0.5 

1 

Case 45: 

0.5 

0 

0.5 

0.5 

Case 46: 

0 

1 

0.5 

0.5 

Case 47: 

0.5 

0 

1 

0.5 

Case 48: 

1 

0.5 

0 

0.5 

Case 49: 

0.5 

0.5 

0.5 

0 

Case 50: 

0.5 

0.5 

0.5 

1 

Case 51: 

1 

1 

1 

0.5 

Case 52: 

0 

0.5 

1 

0.5 

Case 53: 

0.5 

1 

0 

0.5 

Case 54: 

1 

0 

0.5 

0.5 


Table Alb: Testing data. 


TF„,ax = 0.692 + 0.477(0) - 0.687(AHA) - 0.080(AOA) - 0.0650(OPTT) -0.167(0)^ - 0.0129(AHA)(o) + 0.0796(AHA)^ - 
0.0634(AOA)(C^ - 0.0257(AOA)(AHA) + 0.0877(AOA)^ - 0.0521(OPTT)(a) + 0.00156(OPTT)(AHA) + 
0.00198(OPTT)(AOA) + 0.0184(OPTT)l (Al) 

TW 4 = 0.758 + 0.358(0) - 0.807(AHA) + 0.0925(AOA) - 0.0468(OPTT) - 0. 172(o)^ + 0.0106(AHA)(o) + 0.0697(AHA)^ - 
0.146(AOA)(o) - 0.0416(AOA)(AHA) + 0.102(AOA)^ - 0.0694(OPTT)(o) - 0.00503(OPTT)(AHA) + 
0.0151(OPTT)(AOA) + 0.0173(OPTT)l (A2) 

TT„^ = 0.370-0.205(0) + 0.0307(AHA) + 0.108(AOA) + 1.019(OPTT) -O.135(o)^ + 0.0141(AHA)(o) + 0.0998(AHA)^ + 
0.208(AOA)(o) - 0.0301(AOA)(AHA) - 0.226(AOA)^ + 0.353(OPTT)(o) - 0.0497(OPTT)(AOA) - 0.423(OPTT)^ + 
0.202(AHA)(o)^ - 0.281(AOA)(o)^ - 0.342(AHA)^(o) - 0.245 (AHA)^(AOA) + 0.281(AOA)^(AHA) - 0.184(OPTT)\c^ - 
0.281(AHA)( £^(A0A) (A3) 


Xcc = 0.153 - 0.322(0) + 0.396(AHA) + 0.424(AOA) + 0.0226(OPTT) + 0.175(£^^ + 0.0185(AHA)(o) - 0.0701(AHA)^ - 
0.251(AOA)(o) + 0.179(AOA)(AHA) + 0.0150(AOA)^ + 0.0134(OPTT)(c^ + 0.0296(OPTT)(AHA) + 0.0752(OPTT)(AOA) 
+ 0.0192(OPTT)l (A4) 
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